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INTBODUCTION 


Substantial  progress  has  been  made  in  the  last  several  years  in  the  field 
of  statistical  mechanics  of  irreversible  phenomena.  Not  only  the  important 
problem  of  the  irreversibility  Itself  has  been  clarified  to  a  large  extent, 

o 

but  also  it  can  be  claimed  now  that  ve  do  have  a  way  of  describing  kinetic 
equations  describing  particular  physical  situations.3'*5  These  can  be  even* 
tually  solved,  at  least  in  principle,  and  provide  the  transport  coefficients 
of  the  particular  physical  system  under  consideration.  We  refer  in  particular 
to  the  work  of  Prigoglne  and  his  collaborators,  which  has  influenced  the  whole 
field.  '  '  '  Insofar  as  classical  systems  without  external  correlation  sources* 
are  concerned,  the  problem  of  finding  the  kinetic  equations  (for  the  probability 
distributions)  can  be  considered  to  be  solved. ^  On  the  other  hand,  there  Is 
little  doubt  that  the  derivation  of  the  particular  form  of  this  solution,  based 
on  an  infinite  perturbation  expansion  and  the  use  of  Fourier  vectors  in  the 
place  of  the  coordinates,  shall  be  simplified  and  that  alternative  forms  of 
kinetic  equations  still  might  be  found. 

In  view  of  this  generality  of  the  results  of  Prigoglne  and  his  group,  it  ° 
is  Imperative  to  compare  the  other  kinds  of  approaches  with  Prigoglne *s  results; 
such  a  comparison  can  not  only  clarify  some  physical  or  mathematical  aspects  of 
either  theory,  but  also  can  evaluate  the  extent  of  generality  and  the  Inherent 
limitations  involved.  He  attempt  such  a  comparison  with  respect  to  the  theory 
advanced  by  Bogolyoubov  in  19W>. 8 

Bogolyoubov  has  proposed  a  systematic  procedure  of  solving  by  successive 
approximations  the  Bogolyoubov-Born-Oreen-Kirkwood-Yvon  hierarchy  of  integro- 
dlfferential  equations  for  the  probability  distributions  for  sub -sets  of  a 
(s  «N)  particles  of  an  N -particle  system.  Actually  there  are  two  versions  of 
the  Bogolyoubov  theory- -one  uses  the  concentration  C  -  N/v  as  the  expansion 
parameter,  whereas  in  the  other  version  the  potential  energy  of  intexmolecular 
forces,  U,  is  replaced  by  XU  end  an  expansion  in  powers  of  X  is  used.  In  the 
first  version  one  tries  to  reduce  the  problem  of  N-body  motion  to  that  of 
finding  trajectories  and  momenta  predicted  by  an  exact  solution  of  equations  of 
motion  for  2,  3,  ...  s,  bodies.  The  potential  then  has  to  be  assumed  to  be 


*  e.g. ,  turbulence)  the  extent  of  this  limitation  has  not  yet  been  fully 
clarified. 


2 


repulsive.  In  order  to  ensure  that  the  considered  s  bodies  would  be  found  In 
an  uncorrelated  state  (all  relative  distances  very  large)  If  one  goes  to  the 
limit  t— *-«o.  This  version  has  been  considered  by  Choh  and  Uhlenbeck, 
Hollinger,*  Green,*  and  Cohen,*  and  substantial  results  have  been  derived. 

In  the  second  version  the  reduction  to  the  2,  5,  .  .  .  body  problem  Is  ob¬ 
scured  by  the  use  of  a  different  counting  parameter  (l.e. ,  X)  and.  Indeed, 
a  collective  motion  of  an  electron  gas  can  be  described;  on  the  other  hand,  the 
boundary  ccrdltdons  are  much  simpler  and  there  is  no  use  to  assume  that  the 
potential  is  only  repulsive. 

Since  the  Prigogine  theory  is  based  on  a  weak  coupling  expansion,  a  com¬ 
parison  with  the  first  version  would  involve  a  rearrangement  of  Prigogine* s 
equations  so  that  terns  would  be  ordered  with  regard  .to  the  number  of  particles 
Involved  or,  alternatively,  an  expansion  of  Bogolyoubov  results  in  powers  of  X. 
For  that  reason  we  attempted  a  comparison  of  the  weak  coupling  version  of 
Bogolyoubov* s  theory. 

The  weak  coupling  version  has  been  applied  to  the  plasma  problem,  but  It 
has  never  been  formally  solved  to  all  orders  of  X.  Chapter  I  deals  with  this 
problem;  the  Bogolyoubov  procedure  is  strictly  followed  and  the  formal  solution 
Is  given  In  the  form  of  recursion  relations  between  scattering  operators# 

In  a  second  chapter  now  in  preparation  we  express  the  results  of  Chapter  I 
In  terms  of  Fourier  coofeonents  of  the  distribution  function  and  show  how  the 
resulting  expressions  are  related  to  the  kinetic,  equations  of  Prigogine. 


•  References  to  earlier  work  of  this  author  can  also  be  found  In  Chapter  X. 
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CHAPTER  I 


Chapter  I  is  a  reproduction  of  a  paper  by  J .  Stecki,  and  has  been  sub- 
mitted  to  the  Journal  of  Chemical  Physics  for  publication.  _  _ _ 

THE  WEAK  COUPLING  VERSION  OF  BOGOLYOUDOV 1 S  KINETIC  THEORY  OF  GASES*) 

+) 

J  Stecki 

Department  of  Chemistry,  University  of  Southern  California,  Los  Angeles  7, 
Calif. 


Abstract 

The  version  of  Bogolyoubov' s  theory  of  irreversible  phenomena  in  gases, in  which  . 
the  coupling  constant  is  used  as  an  expansion  parameter,  is  generalized  to  all 
orders . 


1 

In  1946  Bogolyoubov  outlined  a  new  approach  to  the  classical  statistical  mecha¬ 
nics  of  irreversible  phenomena  in  gases .  This  theory  has  been  applied  by  Choh  and 
2 

Uhlenbeclc  to  derive  a  generalized  form  of  Boltzmann  equation  in  which  ternary 

3 

collisions  were  taken  into  account.  Recently  Hollinger  has  derived  equations 
valid  to  all  orders  of  concentration, The  generalization  of  Bogolyoubov' s  treet- 

4 

ment  to  all  orders  of  concentration  was  also  considered  by  Cohen  and  Green  and 
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expansion  parameter,  has  also  been  given  by  Bogolyoubov  and  has  been  applied  to  a 

derivation  of  a  kinetic  equation  for  a  homogeneous  electron  plasma  in  a  positive 

6,7,8 

background  .  This  equation  has  also  been  derived  by  Balescu,  significantly,  in 

a  very  different  way,  namely,  by  using  the  theory  developped  by  Prigogine  and  his 
9,10,11  2,4,5,12 
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of  Bogolyoubov  assumptions  and  the  extent  ®f  generality  of  final  kinetic  e^unt^fens 
to  which  his  procedure  leads,  are  important  and  open  questions .  Recently  tftis  pro® 
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blem  has  been  given  preliminary  consideration  bj?5  Re%fhh43  $>  v$o  esmsjdejged  ^rfefe 
2 

ly  the  generalization  of  Boltzmann  equation  in  terms  of  the  recently  npmjUsted’  ge- 
10  1^. 

neral  theory  .  We  intend  to  present  a  discussfofi  §f  this  vefsion  of  Bogfiifyouhov 

®  10 

method  within  the  framework  of  that  general  theory  given  by  Prigogine  .  In  this 
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3.  Stecki  end  H.  S.  Taylor,  to  be  published 


paper  an  explicit  formal  calculation  is  presented  which  leads  to  kinetic  equa¬ 
tions  valid  to  all  orders  of  A-  . 


THE  BOGOLYOGBOV  PROCEDURE  ,f„, 

-v 

Let  us  introduce  the  N-particle  distribution  function  of  momenta,  ,  end 
positions,  q  ,  normalized  to  unity  and  depending  on  time  t.  Let  K  be  the  s-par- 

~i  -  -s  — 

tide  Liouville  operator,  without  the  customary  ^-factors,  defined  as 

Ks 


K°  -  A.  SKs  , 

(1) 

5wK°c;)  > 

(2) 

‘  5jel»l  A#  ' 

(3) 

)*  Vc  --  »  (4>c-s  p;m') 

<*) 

2  u-ij  2L  2iSi/ 

(5) 

%  2p<  o  ?£• 

Oil- 


The  phase  of  the  particle  "i."  is  denoted  by  which  is  understood  to  be  the  6- 

-vector  ^  •  Tliis  is  also  abbreviated  as  (i)  .  Vector  notation  for  p^,  qi? 

v..  is  not  used.  The  operator  K  corresponds  to  a  Poisson  bracket  with  a  hamil- 

— s  ® 

tonian  <§> 

.  .  « —  o.  .  „ 

(S) 


ft-  Z  phi™f  -I-  A  Vi  , 

c  <& 


with  in.temolecu3.ar  potential 


Us  —  zC  U-(j  (|  <},■ ■-  <]jl) 


(?) 


Along  the  mechanical  trajectory  of  the  system,  f^  is  a  constant  of  motion 

+  Kn)  ~fn  ~  (8) 

**) 


and  from  this  the  well-known  B-B-G-K-Y  «  hierarchy  of  equations, 


is  ctaiuvc. ci. 

- 


1,10 


+  Ks)  fa  =  Jdxs+±  ^  ®i,s+i  f<s+j 

C6)i] 


(9) 


Bogolyoubov,  Bora  and  Green,  Kirkwood,  Yvon 


k 


The  reduced  distribution  functions  f  (x, ,  . ..,  x^t)  ere  probability  densities 
for  finding  a  set  of  any  s  particles  at  specified  phases  x^,  . . x^.  The  spatial 

part  of  in  a  homogeneous  system  is  £  =  H/V,  where  V  is  volume  and  £,  number  den- 

s 

sity;  the  spatial  part  of  f^,  reduces  to  £  if  all  particles  are  uncorrelated,  i.  e. 
if  all  relative  distances  between  the  £  particles  are  very  large .  For  this  reason 
it  is  convenient  to  define 

=  c"  (10) 

with  the  restriction  s  <<<  N  insofar  as  the  system  is  finite,  so  that  (9)  takes 


the  form 


(34  +  K?)  fs  =  A  iK5F5  -  fa,  fa-  Bh 


This  system  of  equations  is  not  closed  -  in  contrast  to  macroscopic  transport 
equations .  Bogolyoubov  then  proceeded  to  assume  that  all  higher  distribution  fun¬ 
ctions  are  functionals  of  one-particle  distribution  function,  or  rather  of  the  set 
of  F  '  s  having  positions  and  momenta  of  different  particles  as  arguments;  5^(1, t), 
F  (2,t),  ...  .  Tills  allows  us  to  write 

^  Ts=  0,  \  (12> 

where  0  is  a  time-independent  operator,  and  the  time  £  is  the  same  either  a,s  ar- 

~s 

gument  of  F  '  s  or  of  F  . 

0  — — s  ® 

Suppose  now  one  is  somehow  entitled  to  expand 

f  f]  A  (13) 

a  0 

This  could  be  substituted  into  (ll)  and  the  same  powers  of  A.  collected.  However, 
the  first  of  the  equations  (ll)  reads 

(3t  +  =  Acfae.-tiiM  m 

so  that  to  be  consistent  with  (13)  for  £  =  2  we  have  to  expend  the  time  derivative 


of  Fx  : 


.  f  X*A* 


o 
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with 


AoC0  *  - 


(16) 


A^(0  S  c.  jdx*  0t*  (<■',*)  0*^1).  (17) 

Therefore  we  do  not  expand  F^  itself;  it  follows  that  the  assumed  expansion  (13) 
is  to  he  interpreted  as  the  expansion  of  the  operator  Og  in  powers  of  \  ; 


The  time  derivative  of  F  is  then 


m. 


-X  n?(o  at  _  ?  sf.diV  /w»  - 


a 


SF,(j) 


(17) 


(18) 


or 


»•  ”*»  vM.-«-/vk  <-*  — (rt. 

5tVa  -  X2LX  ,  Us) 

«  O 

where  the  operator  c®m  stands  for  differentiation  with  respect  to  t  followed 
by  substitution  of  _^a(j_)  for  <^.  F^jJ,  for  each  Using  (ll)  and  collecting  the 
same  powers  of  A,  we  arrive  at  the  basic  set  of  equations 

Cp30  +  k°  )  rs°  -  °  > 

_  _wv-i 

(j0  +  k°5)C=  "21  -t-JKsTv'l+  cj<rf*sn£ 


(21) 


(22) 


°  *  on- 1 

which  can  be  solved  successively,  when  following  the  prescription  given  by  Bogo- 

o 

lyoubov.  TOiere  is  a  variety  of  choices  for  the  zeroth  approximation.  Indeed,  the 
equation: 


Ot'  Ki)  f  ”  ° 

has  a  solution 

f  -  fo(h>1r>  -  >  t„)  ,  Cf* *  Cjc  -»;t  J 


(23) 

(24) 


reducing  to  a  prescribed  arbitrary  function  £q(1,  _s;tp)  at  t  =  0.  Guided  by 
the  form  of  the  Boltzmann  equation,  we  want 

%°  *-  IT  VP  •  (25) 

<T' 

This  is  obtained  with  the  aid  of  "boundary  condition"  introduced  by  Bogolyoubov. 

Consider  expression  (24)  which  ca/ also  be  written  as 


-  6  - 

exp  (-«£)  f0 

which  is  noting  else  then  a  Taylor  expansion  of  (2k)  around  t*0  .  The  exponential 
operator  shifts  the  positions  backwards  along  straight  lines  of  unperturbed  notion. 
Then,  for  t  large  enough,  a  situation  will  be  reached,  when  all  the  £  particles  will 
be  distant  from  one  another,  so  that  P  will  reduce  to  a  product  of  P^'s.  So  we  put 


or  in  terms  of 


SU'WI  exp(-rK°s)Fs°  =  TT  ^  tj  V<0  ,  (27) 

dap  (-T  K°s)  T*  *  O  0*^0  .  ^ 

T"*  +  °<=  ' 

These  are  Bogolyoubov  boundary  conditions  for  homogeneous  systems.  In  that  case 
F^  =  aoes  not  depend  on  position.  In  the  inhomogeneous  case  some  care  has 

to  be  exercised,  since  F,  despends  on  a  given  £j  explicitly  end  implicitly  through 
P^i's  which  do  depend  on  positions  now.  The  operator  has  a  property  of  shifting 
everything  which  is  to  its  right,  so  that  F^'s  have  to  be  shifted  separately  in 
the  opposite  direction  in  order  to  compensate  for  the  unwanted  shift  of  their  ar¬ 
guments.  In  Bogolyoubov'  s  notation,  writing  F^x^,  . . .,  x^j  Fp  instead  of  (12), 
we  have 


-fl-vw  &xpC-^K|)Fs  (1,...  =  TT  ^(p.  (2g\ 

T  -»  +  <*>  v  71 


The  exponential  operator  was  denoted  by  Bogolyoubov  by  8^  .  We  split  (29)  into 

^  Ovj*  i£*p(+TK|  ")  ~  ^  1  ^(30) 


iS Ui  <2<p  C-x<s')^rs'C''>”->'SiC!(P('fTl<^T')  “■  0  >  Q"  ^  O.  (31) 

T  -f  00 

These  conditions  have  been  used  by  Bogolyoubov  to  derive  (25)  as  the  solution  of 
the  homogeneous  equation  (21)  and  to  show  that  the  solution  of  an  inhomogeneous 


equation 


-  7  - 


4©o  c*.,- 


(32) 


is 


^  —  ja*  eyLpC-^k*)  C*.  >•••>**>  e*r(+T<W) 


(33) 


whore  g  may  stand  for  tLe  r.h.s.  of  (22). 
s 


FORMAL  S0IOT1CH 


Using  (33)  in  (21)  and  (22)  Bogolyoubov  derived  the  next  approximation,  which 
reads 

}?’  =  \  otr  (-x  *s)  S  '<sfcKt>C^T'cs')  *TT  ^XJ  .  (3h) 

s  }  i  e  y**, 

O 

The  following  features  are  clear.  The  last  term  of  the  r.h.s .  of  (22)  as  well  as 
the  Q0" -operators  through  the  A  -  functions  introduce  l^'s  of  additional  varia¬ 
bles  which  are  always  integrated  over,  as  Indicated  in  (22)  and  in  (17).  To  each 
variable  there  corresponds  a  power  of  c_.  We  denote  these  "dummy  variables"  by  sub¬ 
scripts  <v( ,  o<4  ,  . . .  or  oi  ,  (3  .  The  variables  represent  the  "average  particles 

of  the  medium"  allowed  to  interact  with  the  fixed  set  of  _s  particles  whose  posi- 

n 

tions  and  momenta  are  the  arguments  of  .  For  given  n  there  can  be  at  most  l»n-l 

dummy  variables.  Also  F  's  appear  only  as  a  product.  Therefore  we  assert  that  (l8) 

"1 

can  be  replaced  by 


T«r  -  V*  c'si'JiW  IT  V,(>}  TT  f\  C<)  , 

4«1<5 


ho 


(35) 


_<rv( 

where  the  operators  Als  are  to  be  determined.  It  is  understood  that  they  act 
on  the  main  set  ^  a  ^  and  on  the  set  of  dvmny  variables  <V( ,  . . .  ,  <Vg  .  Hew  the 
results  of  action  of  ^-operators  can  be  written  down  explicitly.  When  making  use 
of  (17)  and  (35),  we  obtain 


-  c0«Trf‘?(O  -  2Z  TT ^Cj)Aa(i) 


u  •  cc-y.)  1  j*\t\ 


(35) 
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i 


p«Pl 


V*0 


(37) 


zy*tfT  K,  (r  *zl  ^F^UimpTr^fs) 


+21  WuriT&on  wmww 

't'O  v',o  55  J  ««W  <*\1)  *  p*Wl  ' 

peWS 

This  term  can  exist  for  n  ^  3.  Finally 


°8,^ 


<n-l 


z  c?  iC  "  Mi  :  21  #  )  otot  it  F.wwTrw. 

£7  5  >  \«t\  (3S) 


We  now  substitute  these  expressions  into  (22)  and.  (33).  We  find  F  in  the  form  (35), 

“8 

which  proves  (35)  by  induction,  (3*0  being  given.  When  collecting  the  some  powers  of 
A  and  £  we  find  three  separate  recursions  -  two  particular  coses 

.xt>(^Kus) 

(39) 


<n,0 

JrO  i^-s 


<*•  ^  ^ 

C'^^s)  S  Ks  _Q.S  nxp  ) 


*='.-1  sT^'  -  <M\,*) -Qzh-.*))  +  (4o) 

-s£'i"x(L*Z  , 

12.”  =  jdr^«p(-xml<?)|SK5P-s 'cxpK^s+c) ^ 

l c ,*))  e  *p(^-  K° + <  ( W )) 

m.'>V’*  O 

,  -w-^v'so  , 

+  Z  )  0v^lVCc,<VxpC^<,((^,o(AVi.KO)j  \  , 

le-SV'l 


and  the  general  one 

I  d  $  Orv-  2. 

<Yl? 


with 


^q'°  =  ^(<uf  c-^>^s)  5Ks«xpO-c,<|)  . 


(42) 


This  is  the  formal  solution,  valid  for  all  n  ^  1.  How  the  operators  SI  can  be  cal¬ 
culated  successively  by  simple  algebra.  Sane  minor  simplificationjls  achieved  by  ma- 
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king  the  integrals  explicit;  there  axe  n  integrations  over  T, ,  ...  ,  xn  end 
t  integrations  over  x  ^ ,  x  p  >  •  •  •  -We  note  that  there  is  a  separate  recursion 
for  l=>  0;  the  operators  J represent  approximations  th  the  s-body  problem  of 
Isolated  s  particles  involved  in  P  .  The  other  separate  recursion  Involves  the  lar- 

**  “TS 

gest  power  of  c  possible  for  given  power  of  A  .  These  terms  have  been  considered 

8 

by  Guernsey  .  For  a  homogeneous  system  they  leed  to  Belescu-Lenard  equation . 

m  HCMXSEKECXJS  SYSTEM 


Considerable  simplifications  occur  if  it  is  assumed  that  does  not  depend  on 
position.  Then 

cxpC-TTK0.)  •  (43) 

Also,  due  to  spherical  symmetry  of  the  interaction  potential 

-  o  .  (44) 

We  can  also  now  identify  the  T- integrals  with  the  free  propagators  familiar  from 

to 

the  scattering  theory,  for  which  many  equivalent  representations  are  known.  Hence 
we  introduce  the  abbreviation 

<*C 

jdx  oxp  (.-tKs  )  =  Q*  .  (45) 

O 

With  (25)  unchanged,  we  obtain 


SI 


,1.0 


9.  s<, 


(46) 


and  the  following  recursions  resulting  from  (39)-(^l) 

•U-V.O 


(-  O 


SI 


m.,o 


&  ks  ns 


m 


4  L  *-*w  {M 

- z  r  nt’X  U 


l  4?  cm.- % 


-  lo¬ 


ll 


<n,rK-  I 


rfx.e(JaSM  (*m-.q»  c^)J 
L 


(49) 


Discussion 

The  last  recursion,  (49),  Is  of  remarkable  simplicity.  It  has  been  examined  by 

8 

Guernsey  In  a  different  but  equivalent  form  and  shewn  to  lead  for  £  a  2  and  when 
used  in  conjunction  with  (15)-(17)  to  the  Baleacu-Lanard  kinetic  equation,  notably 
A25  is  a  sum  of  binary  terms,  so  that  12z(i>Q  Just  cancels  a  term  in  12.^ , 

The  general  term  Is  of  the  form 

§  \*uZ*i\  .... 


so  that  no  two  propagators  follow  one  another  without  being  separated  by  (a  sum  of) 

scattering  operators  6  ^ j .  This  can  be  traced  back  to  the  fact  that  only 

contributes  here.  As  can  be  seen  from  (47),  that  Is  no  longer  true  In  the  general 

term,  'I  4  v  Each  SI  -operator  contains  always  one  ^  at  the  outmost 

left  resulting  directly  from  (33)*  Therefore  the  last  term  of  the  r.h.s.  of  (48) 

will  introduce  a  9,9,  sequence,  which  will  reappear  later  to  give  rise  to  a 

n  o-l 

<3  Q  0tj  sequence,  etc.  We  see  that.  Indeed,  the  terms  with  X  c  are 

exceptionally  simple. 

The  Markovian  character  is  assumed  explicitly  as  early  as  In  (12).  A  non-Marko- 
vien  behaviour  could  result,  if  we  would  allow  the  operator  Og  to  depend  on  time  t^ 
and  it  is  not  clear  hew  this  could  be  done  without  departing  radically  from  the  Bo- 
golyoubov  procedure.  On  the  other  hand  we  know  that  the  most  general  kinetic  equa¬ 
tions  are  non-Markovian1^ . 

Owing  to  the  simplicity  of  the  boundary  conditions,  there  Is  no  need  to  assume 

13  4 

that  the  Intermole cular  potential  is  only  repulsive  ’  *  and  all  expressions  ere  in 
principle  amenable  to  explicit  calculations.  On  the  other  hand,  sane  Individual 
integrals  can  diverge  owing  to  the  repulsive  force  being  infinite  at  the  center  of 


li¬ 


the  particle  for  cannon  potentials;  bmce  sane  res'Jceuvtion  procedure  may  be  needed1^ 
-  a  cannon  feature  of  all  weak- coupling  theories. 
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